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Process and device for determining a transforming element for 
a given transformation function, method and device for 
transforming a digital signal from the time domain into the 
frequency domain and vice versa and computer readable medium 

5 

Cross Reference to Related Applications 

This application claims the benefit of priority of U.S. 
Provisional Application no 60/507,210, filed 29 September 
2003, and U.S. Provisional Application no 60/507,440, filed 
10 29 September 2003, the contents of each being hereby 
incorporated by reference in its entirety for all purposes. 

Further, the following commonly-owned applications are 
concurrently-filed herewith, and herein incorporated in its 
15 entirety: 

"Method for Performing a Domain Transformation of a 
Digital Signal from the Time Domain into the Frequency Domain 
and Vice Versa,'' Atty. Docket No. P100442, and 

"Method for Transforming a Digital Signal form the Time 
20 Domain into the Frequency Domain and Vice Versa," Atty. 
Docket no. P100444. 

This invention relates to a process and a device for 
determining a transforming element for a given transformation 
2.5 f uncti on, a method and a device for transforming a digital 
signal from the time domain into the frequency domain and 
vxce_ srexsa anH a computer readable medium. 

-Domain transformations, for example the discrete cosine 
30 transform (DCT) , are widely used in modern signal processing 
industry. Recently, a variant of the DCT, called integer DCT, 
has attracted a lot of research interests because of its 
important role in' lossless coding applications. The term 
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"lossless" means that the decoder can generate an exact copy 
of the source signal from the encoded bit-stream. 

The DCT is a real-valued block transform. Even if the input 
5 block consists only of integers, the output block of the DCT 
can comprise non-integer components. For convenience, the 
input block is referred to as input vector and the output 
block as output vector. If a vector comprises only integer 
components, it is called an integer vector. In contrast to 

10 the DCT, the integer DCT generates an integer output vector 
from an integer input vector. For the same integer input 
vector, the integer output vector of integer DCT closely 
approximates the real output vector of DCT. Thus the integer 
DCT keeps all the good properties of the DCT in spectrum 

15 analysis. 

An important property of the integer DCT is reversibility. 

Reversibility means that there exists an integer inverse DCT 

. 

(IDCT) so that if the integer DCT generates an output vector 

20 £ from an input vector x, the integer IDCT can recover the 

vector x from the vector jf. Sometimes the integer DCT is also 

referred to as the forward transform, and the integer IDCT as 
the backward or inverse transform. 

25 A transform called integer modified discrete cosine transform 
(IntMDCT) is recently proposed and used in the ISO/IEC MPEG-4 
audio compression. The IntMDCT can be derived from its 
prototype - the modified discrete cosine transform (MDCT) . In 
[1], Malvar gave an efficient realization of MDCT by 

30 cascading a bank of Givens rotations with a DCT-IV block. It 
is well known that Givens rotation can be factorised into 
three lifting steps for mapping integers to integers, see for 
example [2] . 
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Therefore, the realization of IntMDCT relies on an efficient 
implementation of integer DCT-IV. 

5 Integer transforms can be directly converted from their 
prototypes by replacing each Givens rotation with three 
lifting steps. Because in each lifting step there is one 
rounding operation, the total rounding number of an integer 
transform is three times the Givens rotation number of the 
10 prototype transform. For discrete trigonometric transforms 
(for example the Discrete Fourier Transform (DFT) or the 
Discrete Cosine Transform (DCT) ) , the number of Givens 
rotations involved is usual.ly at Nlog 2 N level, where N is 

the size of the blocks , i.e. the amount of data symbols 
15 included in each block, the digital signal is divided into. 

Accordingly, the total rounding number is also at Nlog 2 N 

level for the family of directly converted integer 
transforms. Because of the roundings, an integer transform 
only approximates its floating-point prototype. The 
20 approximation error increases with the number of rounding 
operations . 

The invention solves the problem of determining a 
transforming element for a given transformation function, 

25 which transformation function comprises a transformation 
matrix and corresponds to a transformation of a digital 
signal from the time domain into the . frequency domain or vice 
VtnrsHT — such that the number of roundings comprised by the 
transforming element is significantly reduced. The invention 

30 further provides a method for transforming a digital signal 
f"roTu~the~ time domain into the frequency domain or vice versa 
according to the determined transforming element. 
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The problem is solved by a process and a device for 
determining a transforming element for a given transformation 
function, a method and a device for transforming a digital 
signal from the time domain into the frequency domain or vice 
5 versa and a computer readable medium with the features 
according to the independent claims. 

According to the invention, there is provided a process for 
determining a transforming element for a given transformation 

10 function, which transformation function comprises a 

transformation matrix and corresponds to a transformation of 
a digital signal from the time domain into the frequency 
domain or vice versa, wherein the transformation matrix is 
decomposed into a rotation matrix and an auxiliary matrix 

15 which, when multiplied with itself, equals a permutation 

r 

matrix multiplied with an integer diagonal matrix; the 

■ 

rotation matrix and the auxiliary matrix are each decomposed 
into a plurality of lifting matrices; and the transforming 
element is determined to comprise of a plurality of lifting 
20 stages which correspond to the lifting matrices. 

Further, according to the invention there is provided a 
device adapted to perform the process described above. 

25 Further, according to the invention, there is provided a 

method for transforming a digital signal from the time domain 
into the frequency domain or vice versa using a transforming 
element, wherein the transforming element corresponds to a 
given transformation function, which transformation function . 

30 comprises a transformation matrix wherein the transforming 

element is determined by a process comprising decomposing the 
transformation matrix into a rotation matrix and an auxiliary 
matrix which, when multiplied with itself, equals a 
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permutation matrix multiplied with an integer diagonal 
matrix; decomposing the rotation matrix and the auxiliary 
matrix each into a plurality of lifting matrices; and 
determining the transforming element to comprise of a 
5 plurality of lifting stages which correspond to the lifting 
matrices; and wherein further each lifting stage comprises 
the processing of sub-blocks of the digital signal by an 
auxiliary transformation and by a rounding unit. 

10 Further, according to the invention there is provided a 
device adapted to perform the method described above. 

Further, according to the invention, there is provided a 
computer readable medium having a program recorded thereon, 
wherein the program is adapted to make a computer perform a 
process for determining a transforming element for a given 
transformation function, which transformation function 
comprises a transformation matrix and corresponds to a 
transformation of a digital signal from the time domain into 
the frequency domain or vice versa, wherein the 
transformation matrix is decomposed into a rotation matrix 
and an auxiliary matrix which, when multiplied with itself, 
equals a permutation matrix multiplied with an integer 
diagonal matrix; the rotation matrix and the auxiliary matrix 
are each decomposed into a plurality of lifting matrices; and 

* 

the transforming element is determined to comprise of a 
plurality of lifting stages which correspond to the lifting 
matrices . 

30 Further, according to the invention, there is provided a 

computer readable medium having a program recorded thereon, 
wherein the program is adapted to make a computer perform a 
meth od for transformi ng a digital signal from the time domain 
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into the frequency domain or vice versa using a transforming 
element, wherein the transforming element corresponds to a 
given transformation function, which transformation function 
comprises a transformation matrix wherein the transforming 
5 element is determined by a process comprising decomposing the 
transformation matrix into a rotation matrix and an auxiliary 
matrix which, when multiplied with itself, equals a 
permutation matrix multiplied with an integer diagonal 
matrix; decomposing the rotation matrix and the auxiliary 
10 matrix each into a plurality of lifting matrices; determining 
the transforming element to comprise of a plurality of 
lifting stages which correspond to the lifting matrices; and 
wherein further each lifting stage comprises the processing 

* 

of sub-blocks of the digital signal by auxiliary 
15 transformations and by a rounding unit. 

In a preferred embodiment, the invention provides a process 
and a method for realizing the integer type-IV DCT 
transformation. The method according to the invention 
requires a significantly reduced number of roundings compared 
to prior art methods. As a result, the approximation error is 
greatly reduced, in the case of DCT-IV it can be reduced from 
the usual Nlog 2 N level to be as low as 2.52V, where N denotes 
the block size of the digital signal. The method according to 
the invention is low in computational complexity and modular 
in structure . 

The method and the device according to the invention can be 
used for any types of digital signals, such as audio, image 
30 or video signals. The digital signal, which is a digitised 
signal, corresponds to a physical measured signal, may be 
generated by scanning at least one characteristic feature of 
a corresponding analog signal (for example, the luminance and 
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chrominance values of a video signal , the amplitude of an 
analog sound signal, or the analog sensing signal from a 
sensor) . The digital signal comprises a plurality of data 
symbols. The data symbols of the digital signal are grouped 
into blocks , with each block having the same predefined 
number of data symbols based on the sampling rate of the 
corresponding analog signal. 

The method according to the invention can be used for 
transforming an input digital signal which represents integer 
values to an output signal which also represents integer 
values. The transformation method according to the invention 
is reversible. The output signal may be transformed back to 
the original input signal by performing the transformation 
method according to the invention. Such a reversibility 
property of the transformation according to the method of the 
invention can be used in lossless coding in which the output 
signal should be identical to the original input signal. 

Such an integer transformation of signals according to the 
invention can be used in many applications and systems such 
as MPEG audio, image and video compression, JPEG2000 or 
spectral analyzers (for analyzing Infrared, Ultra-violet or 
Nuclear Magnetic Radiation signals) . It can also be easily 
implemented in hardware systems such as in a fixed-point 
Digital Signal Processor (DSP), without having to consider 
factors such as overflow in the case of a real— value signal 
transformation . 

According to the method according to the invention, the 
digital signal is transformed to the frequency domain by a 
transforming element, which is according to the process 
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according to the invention determined for a given 
transformation function. 

Preferably, the transforming element comprises a plurality of 
lifting stages. 

The transforming element can be illustrated based on the 
model of a lifting ladder. The lifting ladder model has two 
side pieces, each for receiving one of two groups of data 
symbols. Two or more cascading lifting stages are provided 
between the two side pieces. Each lifting stage receives a 
signal at one end (input end), and outputs a signal at the 
other end (output end) via a summation unit. A rounding unit 
is arranged at the output end. The lifting stages are 
arranged between the side pieces in an alternating manner, 
such that the output (or input) ends of adjacent lifting 
stages are connected to the different side pieces . 

It should be noted that although the transforming element is 
described in the form of the lifting ladder model, it is only 
to illustrate the transformation paths of the transforming 
element. However, the invention shall not be limited to said 
ladder model . 

The number of lifting stages of the transforming element is 
defined by the number of lifting matrices which is determined 
by the process according to the invention. 

Discrete cosine transforms, discrete sine transforms, 
discrete Fourier transforms or discrete W transforms are 
examples of transformation functions that may be used as the 
transformation function according to the invention. The 
number of lifting sLayes -af- the tran sf orming element may be 
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different, depending on the result of the process according 
to the invention for determining a transforming element for 
the respective transformation functions. 

5 Illustrative embodiments of the invention are explained below 
with reference to the drawings, in which 

Figure 1 shows the architecture of an audio encoder according 
to an embodiment of the invention. 

10 

Figure 2 shows the architecture of an audio decoder according 
to an embodiment of the invention, which corresponds to the 
audio coder shown in figure 1. 

Figure 3 illustrates an embodiment of the process according 
to the invention, wherein the transformation function is a 
DCT-IV transformation function. 

■ 

Figure 4 shows a flow chart of an embodiment of the method 
according to the invention. 

Figure 5 illustrates an embodiment of the method according to 
the invention using DCT-IV as the transformation function. 

» 

Figure 6 illustrates the algorithm for the reverse 
transformation according to the embodiment of the method 
according to the invention illustrated in figure 5. 

Figure 7 shows the architecture for an image archiving system 
30 according to an embodiment of the invention. 

Elaure 8 illustrates an embodiment of the method according to 
the invention using- DWT-IV crs—tire tr ansforma tion function. 
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Figure 9 illustrates the algorithm for the reverse 
transformation according to the embodiment of the method 
according to the invention illustrated in figure 8. 

* 

5 

Figure 1 shows the architecture of an audio encoder 100 
according to an embodiment of the invention. 

The audio encoder 100 comprises a conventional perceptual 
10 base layer coder based on the modified discrete cosine 

transform (MDCT) and a lossless enhancement coder based on 
the integer modified discrete cosine transform (IntMDCT) . 

An audio signal 109 which, for instance, is provided by a 
15 microphone 110 and which is digitalized by a analog-to- 
digital converter 111 is provided to the audio encoder 100* 
The audio signal 109 comprises a plurality of data symbols. 
The audio signal 109 is divided into a plurality of blocks, 
wherein each block comprises a plurality of data symbols of 
20 the digital signal and each block is transformed by a 

modified discrete cosine transform (MDCT) device 101. The 
MDCT coefficients provided by the MDCT device 101 are 
quantized by a quantizer 103 with the help of a perceptual 
model 102. The perceptual model controls the quantizer 103 in 
25 such a way that the audible distortions resulting from the 
quantization error are low. The output of the quantizer 103 
is encoded by a bitstream encoder 104 which produces the 
lossy perceptually coded output bitstream 112. The bitstream 
encoder 104 losslessly compresses its input to produce an 
30 output which has a lower average bit-rate than its input by 

standard methods such as Huffman-Coding or Run-Length-Coding. 
The input audio signal 109 is also fed into an IntMDCT device 
105 which produces IntMDCT coe fficients . The quantized MDCT 
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coefficients, which are the output of the quantizer 103, are 
used to predict the IntMDCT coefficients. The quantized MDCT 
coefficients are fed into an inverse quantizer 106 and the 
output of the inverse quantizer 106 is fed into a rounding 
5 unit 107 by which they are rounded to integer values and the 
residual IntMDCT coefficients, which are the difference 
between the output of the rounding unit 107 and the IntMDCT 
coefficients, are entropy coded by an entropy coder 108, 
which, analogous to the bitstream encoder 104, losslessly 
10 reduces the average bit-rate of its input and produces a 

lossless enhancement bitstream 113, which, together with the 
perceptually coded bitstream 112, carries the necessary 
information to reconstruct the input audio signal 109 
exactly. 

15 

Figure 2 shows the architecture of an audio decoder 200 
comprising an embodiment of the invention, which corresponds 
to the audio coder 100 shown in figure 1. 



20 The perceptually coded bitstream 207 is decoded by a 

bitstream decoder 201, which performs the inverse operations 
to the operations of the bitstream encoder 104 of figure 1 
and is fed to an inverse quantizer 202. To its output, the 
inverse MDCT is applied by an inverse MDCT device 203. Thus, 

25 the reconstructed perceptually coded audio signal 209 is 

obtained. The lossless enhancement bitstream 208 is decoded 
by an entropy decoder 204, which performs the inverse 
operations to the operations of the entropy encoder 108 of 
figure 1 and which produces the corresponding residual 

30 IntMDCT coefficients. The output of the inverse quantizer 202 
is rounded by a rounding device 205 and is added to the 
residual IntMDCT coefficients, thus producing the IntMDCT 
coefficients. Finally, the inverse IntMDCT is applied to thre 
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IntMDCT coefficients by an inverse IntMDCT device 206 to 
produce the reconstructed losslessly coded audio signal 210. 



As mentioned above, it is shown in [2] that the core of 
5 IntMDCT, which plays an important role in lossless audio 

coding and which is used in the embodiment of the invention 
illustrated in figures 1 and 2, is an integer DCT-IV. 

Figure 3 illustrates a flowchart of an embodiment of the 
10 process according to the invention, wherein the 

transformation function is a DCT-IV transformation function. 

In the following, an embodiment of the process according to 
the invention for determining a transforming element for a 

15 DCT-IV transformation function is explained. The determined 

transforming element is used in the encoder shown in figure 1 
to implement IntMDCT and the according inverse transforming 
element is used in the decoder shown in figure 2 to implement 
inverse IntMDCT. For a description how to implement IntMDCT 

20 and inverse IntMDCT with DCT-IV, see [2] . 

The DCT-IV transformation function with a tf-point real input 
sequence x(n) is defined as follows (see [2]): 




(1) 



Let C™ be. the transformation matrix of DCT-IV, that is 




(2) 



30 
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According to the embodiment of process according to the 
invention, the transformation matrix C™ is decomposed into 

a rotation matrix and a matrix which, when multiplied- with 
itself, equals a permutation matrix multiplied with an 
integer diagonal matrix. 



For clarity, N is assumed to be even 



In step 300, the process is started. 



In step 301 the terms of the transformation function for even 
indices are separated from the terms with odd indices: 



HT / 1 1 n \ 



^W»)«x>s(0n + f)C2» + f)i) 

/ 



I ~-l 

^ ^ x(2n + 1) cos((m + |)(2n + 1 + -) i j 



^-^(„)cos 



(f) 



r— 1-1 

J — y x 2 (n)cos 



(/« + — )(«+—)— ■ 



(f) 



thus, 



y(m) = 



jCj(n)cos 



(m + -)(n + -) 



(f) 



— (m + — ) 



n 



(3) 
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N 

where xi with the components x 1 (n) = x(2n) ,/i = 0,...,— -1 is a 
vector comprising all x (n) with even indices and x 2 with the 



N 

components x 2 (n) - x(2n + 1) ,n = 0 v ..,— -1 is a vector comprising 
all x(n) with odd indices. 



The following two abbreviations are used: 



a mn =(m + — + — ,m,#i-0,...JV-l (4) 

2 2 



(f) 



^ m= I(m + i)-^ 7 = (m-fi)^ 7 ,m = (W\T-l (5) 
w 4 2 / N\ 2 2N 



With equations (4) and (5), equation (3) can be written as 



yim^^^^cosia^-pj + gx 2 (/i)cos(a^+^ m ) (6) 

In step 302, the following two addition theorems for the 
cosine are applied: 

cos(a ro ^ ~P m ) = cosa^ ■ cos p m + sin -san/^ ( 7 ) 

cos(a m ^ + £ J = cosa^ -cos/3^ -sina w>/1 -sin^ (8) 



With (7) and (8), equation (6) can be written as 
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2 



2 



N 





Jrj^^C^cos^ +^ 2 (/i)cosj3 nj ]cosa JD>n + ^[^(^sin^ - x 2 (n) sin Jsin a 




( 



2 



N 




^^ SL [x 1 (n)+x 2 (n)]cosa m>n + g^^k(n)-x 2 (n)]sina mi) 



(9) 



5 for abbreviation, the two N/2-dimensional vectors x 2 ' and xs' 
are defined to comprise the following components: 



^\n) = -j=[x 1 (n) + x 2 (n)] ,n~0,...,— -l 

1 

x 2 '(«) - j^k (») - ^2 (»)] > w - y " 1 



10 



With (10) and . (11), equation (9) simplifies to 



(10) 



(11) 



rtr 



y(m) = 



(t) 



TV 




1 

cos^Xj'C^cosa^ + Y sin 0„V(*) sin 



(12) 



15 In step 303, the vector 



ry(0) 



is divided into the two parts y.o anc * Yi wherein 

20 
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y 0 



y0) 



(13) 



yi 



y(N-l) 
y(N-2) 



y(-) 



(14) 



5 The vector vj comprises the components of j£ corresponding to 

N 

the indices from — to N -1 in reverse order. 

2 



The components y x (ni) of the vector fulfil the following 
equation: 



10 



y^m) = y(N -1- m) 



N 



ra,n 



(15) 



for in =»0,...,— — 1 . 

£0 



15 Note that 



N 



(f ) I s 



2_ 2 

cos /J ^ JCj r (n) cos a N + V sin p N x 2 r (/i) sin a N 



2 2 " 0 2 



(16) 



for — <UU. 
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In step 304, %p and yi are each expressed by a DCT-IV matrix 

N 

C% and a DST-IV matrix each of dimension — . 



5 This is achieved in the following way: 



1 K 

For B =(m + — ) the following equations hold 

Hm v 2 J 2N 



B w =(m + - + — ) = (m + — ) — - + — ■-— = j3 m +— (17) 

«4 V 2 2 2N v 2'2N 2 2N 4 



Further, 



cos =cos(|-/U = sin (19) 
15 sin =sin(|-/U = cos 0 m < 20 > 



1 1 7T 

For a mtn = (m + — )(/* + — ) the following equation holds: 

2 2 



(f) 



20 a 



1. 7T 



N-i-m.il v 2 /v 2 (N\ 2 2 / J\T 



(f) 



= {N -m- -)(« + -) 



1. 



1 11 1 

= ^(n+-)y^-(m + -)(n + -)y^- = 2w(n + -)-a ))Ijn = 2*r n + * - 



(f) 



(21) 
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Further: 



cosa^ m „ - cos(2tt n + n-a m J = cos(tt -a m J = -cos a 



sina N _. = sin(27r n + it - ) = sin(?r - a^) = sina ra> , 



(22) 
(23) 



With the equations (19), (20), (22) and (23) equation (15) 
can be rewritten: 



y(N-l-m) = 



(f) 



'(^cosa^.j.^ + ^ sin P N ^ m x 2 1 (n) sin a N ^_ mtn 

n 



n 



CN 



10 



ll 



(f) 



AT 



^sin )5 m jc, '(«)(- cos a miB ) + £ cos /3 m * 2 '(^sina^ 

n n 



N 



N 



i 



= -sin£ m 



1 



cos a 



+ cos0 m 



(n) sin a 



/re ,n 



(24) 



for m = 0,...,— -1 - 

2 



15 



Since further equation (12) yields 



s 



y(m) 



i 



J V cos j3 m *! '(n) cosa mn + Y sin fl. x 2 '(«) sin a mn 
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= <*> S An 



N 



cos a 



N 



(fr 



(25) 



the expressions for £o and £i can be formed from (24) and 
(25): 



y 0 - ttiag[cos p m } C% ■ ' + gggjsin fl» ] • • x 2 ' , m=0,...,^-l 

2 2 



(26) 



=£//ag[-sin/3 m }c^ • x/ + diag[cos p m ]• 5 



2 



n ^ 1 

/W = U,..., JL 

2 



(27) 



■ • 

where denotes a N/2-by-N/2 diagonal matrix with a m in 

10 the mth row, C# is the transformation matrix of the DCT-IV 



transformation and S„ is the transformation matrix of the 



N_ 
2 



discrete sine transformation of type four (DST-IV) 



The- two equations (26) and (27) can be expressed as a single 
15 one in the following way: 



y 0 (0) 



.v fft\ 



( 1S 1\ 



sin fi 1 



—sin p x 



cos >3 2 



-•sin/? m 



>IV 

2 



2 



V(f-i) 
V(0) 



(28) 



N 

20 In the following, the — x — - matrices 

2 2 
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J„ - 



N 

2 



(29) 



and 



N_ 
2 



(30) 



are used. 



10 One sees that 



y(2) 

* 

y(N-l) 



N 
2 



N_ 
2 



yo(r 2 ^ 
*(0) 



(31) 



which is abbreviated by 



15 



2 



2 



(32) 



One further sees that with the Nx-N -aiaLiix defi ned as 
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In/2 In/2 
I -I 

J N/2 J N/2 



(33) 



the equation 



" V(0) • 

• 






• 

• 

V(f-D 
V(0) 

* 




• 

*(f-l) 

x 2 (0) 

• 


• 




• 



(34) 



holds 



10 Equation (34) is abbreviated by 



R 



(35) 



Further, let be an even-odd matrix, i.e., a permutation 
15 matrix reordering the components of the vector 



x = 



*(0) 



x(N -1) 



by separa ting the components which correspond to even indices 
2.0- frnm the components which correspond to odd indices, i.e., 



•lit 



•111 



WO 2005/031597 



PCT/SG2004/000122 



22 



* 2 (0) 



f N f\ 



= p 



eo 



x(0) 



x(N-l) 



(36) 



or shortly 



= P„x 



(37) 



holds . 



By combining the equation (28) with the equations (31) , (34) 
10 and (36) one gets: 



y(2) 



y(N-l) 



2 



N 

2 



cosfi t 



COS/^ 



2 



*(0) 



15 



with the above abbreviations and the abbreviations 



(38) 



"cos/?! 



C = 



(39) 



20 and 
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sin ft 



(40) 



equation (38) can be shortly written as 



N 
2 



2 



c 



pr eg — 



(41) 



IV 



In step 305 

1 



is expressed by . This is done with the 



equation: 



= 7 -C™ D 

2 2 2 2 



(42) 



wherein D Nn is a diagonal matrix of dimension N/2 given by 



-1 



D„ = 



N 
2 



(43) 



-1 



With equation (42) equation (41) can be written as 
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In step 306 the NxN rotation matrix is calculated, 

which comprises the first three matrices of equation (44) 



PO 



2 



N_ 
2 



C SI 
-5 C 



2 



2 



(45) 



Multiplication of the three matrices in equation (45) yields 



R po 5=8 



It 



COS 



AN 



cos 



3tt 
AN 



— sm 



3n 
AN 



— sin 



4iV 



cos 



SIB 



AN 

. (AT-l>r 
-sin 

AN 



cos 



4JV 
4tf 



it 



sin 



4tf 



sin 



3?r 
AN 



cos 



3tt 
4JV 



7r 



cos 



AN 



(46) 



In step 307, the auxiliary matrix, which, when multiplied 
with itself, equals a permutation matrix multiplied with an 
integer diagonal matrix, is calculated. 

The auxiliary matrix comprises the fourth and fifth matrices 
of equation (44) : 



^N/Z 



C™ D 



J2 



C™ D 



-C™ D 

^niz-^niz 



(47) 



Note that 
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1 0 
0 1 



0 

0 1 
-1 0 



0 



-1 



25 



0 



0 



2 



D 



2 



0 



0 



0 



0 



0 



0 



0 



0 



(48) 



I.e., T multiplied with itself equals a permutation matrix 
5 multiplied with an integer diagonal matrix. 



With equations (46) and (47), equation (44) simplifies to 



10 



(49) 



Thus, the transformation matrix C™ is decomposed into a 

rotation matrix , into an auxiliary matrix T that, when 

multiplied with itself, equals a permutation matrix 
multiplied with an integer diagonal matrix and into an even 



15 odd matrix R, 



eo 



In step 308 and T are each factorised into a product of 



lifting matrices. T is factorised in the following way: 



20 r = j a r 2 r 3 = 



1 N/2 



N/2 



-z> 



N/2 



N/2 



N/2 



N/2 



(50) 



where 
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(51) 



K 3 = 'yf^C^ /2 D N/2 + In/2 



(52) 



(53) 



And R is factorised according to the following equation: 
p° 



R po = RiR 2 R 3 = 



JV/2 



JV/2 



#/2 f*2 



JV/2 



*N/2 



(54) 



where 



H,=H, 



- tan 



3;r 



- tan 



SN 



- tan 



(N - l);r 
8iV 



(55) 



and 



H 2 = 



n 



sin 



4JV 



sin 



3tt 
4JV 



sin 

4JV 



(56) 



Thus, Equation (4 9) can be re-written as 



C# - R\R2^Z^\^2^Z^ id 



(57) 
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In step 309, the lifting matrices are aggregated as far as 
possible. Matrix S is defined as the product of and T\, 

that is 



5 5 =^3^ - 



N/2 



N/2 



(58) 



This aggregation is possible since S is also a lifting 
matrix. 

10 From (57) and (58), the final factorisation formula for DCT- 
IV matrix is obtained as 



C% =^ 2 5 r 2 r 3 P go (59) 



15 Equation (59) indicates that the transforming element for the 
integer DCT-IV transformation according to the invention 
comprises five lifting stages. 



Since the final factorisation formula is* determined, the 
20 process is stopped in step 310. 

Figure 4 shows a flow chart 400 of an embodiment of the 
method according to the invention using five lifting stages, 
a- Xlnst lifting stage 401, a second lifting stage 402, a . 

25 third lifting stage 403, a fourth lifting stage 404 and a 
fifth lifting stage 4 05 . tm « met h od is prefexably used in 
the IrrtMDCT device 105 of figure 1 and the inverse IntMDCT 
device^ 2JD£ of figure 2 to implement IntMDCT and inverse 
IntMDCT ,- respectively. -In. .Fig— 4., .andL.xg .axe first and 

30 second blocks of the digital signal, respectively. z lf z 2 , z 3 
are intermediate signals, and yi and are output signals 
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corresponding to the first and second block of the digital 
signal, respectively. 

Figure 5 shows a flow chart of an embodiment of the method 
5 according to the invention, wherein the transformation 
function is a DCT-IV transformation function. The 
transforming element used in this embodiment corresponds to 
the equation (59), i.e., is the one which is determined by 
the embodiment of the process illustrated in figure 3. 

10 

The transforming element comprises five lifting stages which 
correspond to the five lifting matrices of equation (59) . 

Further, the transforming element comprises a data shuffling 
15 stage corresponding to the permutation matrix . 

In figure 5, the input of the first lifting stage are the two 
blocks of the digital signal xi anc * i£2, iEi/ *z anci are 
intermediate signals and and y.2 are output signals 
20 corresponding to the first and second block of the digital 
signal , respectively . 

The input of the transforming element x and the two input 
blocks of the first lifting stage of the transforming element 
25 Ki and x 2 fulfil the equation 




(60) 



(confer equation (37)) 

30 
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In the following, the first lifting stage 501 is explained, 
which is the lifting stage corresponding to the lifting 

matrix T, . 



Let 



be the output vector of the first lifting stage, for 



the time being without rounding to integer values, i.e 



(61) 



Using the definition of 7^ provided by equation (50), 
equation (61) is re-written as 



JV 
2 



jST, I 



2 



(62) 



Since in this embodiment, a reversible algorithm for integer 
DCT-IV in this provided, there is included a rounding to 
integer values. So, according to equation (62), in the first 
step 506 of the first lifting stage 501, xj is multiplied by 
. The result of this multiplication is rounded to integer 

values in step 507. The rounded values are then added to x 2 
in step 508. Thus, the intermediate signal zi fulfils the 
equation 



(63) 



where [*] denotes rounding operation 
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Since the other four lifting stages 502, 503, 504, 505 of the 

transforming element illustrated in figure 5, which 

correspond to the matrices 7^ , S, and respectively, 

have identical structures as the first lifting stage 501, 
5 their explanation is omitted. It only has to be noted that in 
the adding step 509 of the second lifting step 502, xj is, 
according to the definition of 7^ , multiplied with - D N/2 • 

In the following, the lifting stages of the transforming 
10 element of the reverse transformation are explained with 
reference to figure 6. 

Figure 6 illustrates the lifting stages of the transforming 
element for the reverse transformation of the transformation 
15 illustrated in figure 5. 

In figure 6, the input of the first lifting stage are the two 
blocks of the digital signal j£j and j/ 2 , Zx, z 2 and z 3 are 
intermediate signals and x 2 and x 2 are output signals 
20 corresponding to the first and second block of the digital 
signal, respectively. 

The last lifting stage 605 illustrated in figure 6 is inverse 
to the first lifting stage 501 illustrated in figure 5. So, 
25 in the first step 606 of the last lifting stage 605, x 2 is 
multiplied by K 3 . The result of this multiplication is 

rounded to integer values in step 607. The rounded values are 
then subtracted from z 2 in step 608. Thus, the signal x 2 
fulfils the equation 



30 
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where |*J denotes rounding operation. 

Since the other four lifting stages 601, 602, 603, 604 of the 
transforming element illustrated in figure 6, which are 
inverse to the lifting stages 505, 504, 503 and 502 
respectively, have identical structures as the last lifting 
stage 605, their explanation is omitted. It only has to be 
noted that after the adding step 609 of the fourth lifting 
step 604, the result of the adding step 609 is multiplied 
with -D Nn to produce Xj. 

It can be seen that the lifting stages 605, 604, 603, 602 and 
601 of figure 6 are inverse to the lifting stages 501 to 505 
of figure 5, respectively. Since also the. permutation of the 
input signal corresponding to the matrix can be reversed 

and an according data shuffling stage is comprised by the 
inverse transforming element, the provided method is 
reversible, thus, if used in the audio encoder 100 and the 
audio decoder 200 illustrated in figures 1 and 2, providing a 
method and an apparatus for lossless audio coding. 

An analysis of the number of roundings used in this 
embodiment is given at the end of the description of the 
..invention. 

Figure 7 shows the architecture for an image archiving system 
acooxdxnrg to an embodiment of the invention. 

In figure 7 an image source 701, for instance a camera, 
provides an analog image signal. The image signal is 
processed by a analog- to-digital converter 702 to provide a 
corresponding digital image signal. The digital image signal 
is losslessy encoded by a lossless image encoder 703 which 
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includes a transformation from the time domain to the 
frequency domain. In this embodiment, the time domain 
corresponds to the coordinate space of the image. The 
lossless coded image signal is stored in a storage device 
704, for example a hard disk or a DVD. When the image is 
needed, the losslessly coded image signal is fetched from the 
storage device 704 and provided to a lossless image decoder 
705 which decodes the losslessly coded image signal and 
reconstructs the original image signal without any data loss. 

Such lossless archiving of image signals is important, for 
example, in the case that the images are error maps of 
semiconductor wafers and have to be stored for later 
analysis . 

In the following, a further embodiment of the method for the 
transformation of a digital signal from the time domain to 
the frequency domain and vice versa according to the 
invention is explained. This embodiment is preferably used in 
the lossless image encoder 703 and the lossless image decoder 
705 of the image archiving system illustrated in figure 7. 

Figure 8 illustrates an embodiment of the method according to 
the invention using DWT-IV as the transformation function. 

The DWT-IV of a N-point real input sequence x(ri) is defined 
as follows: 

y(m) - ^— 2*(*)sinf — + N — J m > n = ^* " iN ( * 

The transformation matrix of DWT-IV is given by 
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. / n (m + y2)(n+l/2)27C 

sin — + — — — - — 

^4 N 




(66) 



04,— ,N-1 



According to the process according to the invention, the DWT 
IV matrix is factorised into the following form: 



«N V = RnT % 



(67) 



Rjf is a NxN rotation matrix defined as 



1 0 iv v = — f= 
— V2 



N/2 J N/2 



~~Jn/2 In/2 



(68) 



15 



I Nf2 is the identity matrix of order N/2 (confer equation 
(29) ) • J N/2 is the counter identity matrix of order N/2 
(confer equation (30)). 



P N is a JVxJV permutation matrix 



N 



I 



N/2 



N/2 



(69) 



20 T is a JVxJVmatrix given by 



, ""N/2 

sf2 Cn/2&N/2 



D 

K -N/2 xy N/2 



(70) 



where C" n is the DCT-IV matrix of order ST/2 

25 
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^N/2 



(N/2) 



cos 



(m+]/2)(n +l/2)n 
(N/2) 



(71) 



D N/2 is an order N/2 diagonal matrix given by 



5 D 



N/2 



-1 



-1 



(72) 



According to the process according to the invention, R^ and 
T are further factorised into a product of lifting matrices: 



10 T = T t T 2 T 3 



N/2 



N/2 



D 



N/2 **"2 
J 



N/2 



\I 



N/2 



^3 In/2 



(73) 



where K t 



= (V2 



^n/2 "^Nn^Nn^P 



N/2 *^2 



^Nf 



V2 



K 3 — -yf2C# /2 D N/2 J N/2 , 



and 



15 R N =*R t R 2 R 3 = 



AT/2 



Jtf/2 



AT/2 



'N/2 

"3 In/2 



(74) 



where 



=H 3 =-tan(n/8)-J 



N/2 



- 0.414 



- 0.414 



-0.414 



(75) 



20 



and 
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0.707 



0.707 



0.707 



(76) 



Thus, equation (67) can be written in the form 



W^ V = RiR 2 R3TiT 2 T3% 



(77) 



According to the process according to the invention, the 
lifting steps are aggregated as far as possible. 



10 In this embodiment the lifting matrices R 3 and T t can be 
aggregated to the lifting matrix 5 : 



S = R 3 T t 



N/2 



N/2 



(78) 



15 From (77). and (78) , the following factorisation formula for 
the DWT-IV matrix is obtained: 



Wn V - RlR2§ T2T3% 



(79) 



20 Equation (79) indicates that the transforming element for the 
integer DWT-IV transform according to the invention comprises 
five lifting stages. 



Further, the transforming element comprises a data shuffling 
25 stage corresponding to the permutation matrix P N . The data 

shuffling stage rearranges the components order in each input 
data block. According to P N , the input data vector is 

rearranged in the following way: the first half of vector 
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remains unchanged; the second half of vector is flipped, 



i.e. 



*1 




*1 


• 
• 




• 


* 

X Nt2 




• 

X N/2 


X N/2+l 




X N 






X N-1 


X N-1 




» 
* 

• 






X N/2+l 



(80) 



10 



In figure 8, the input of the first lifting stage are the two 
blocks of the digital signal xi and *2, zi* z 2 and z 3 are 
intermediate signals and j/i and % 2 are output signals 
corresponding to the first and second block of the digital 
signal, respectively. 



15 



The input of the transforming element x and the two input 
blocks of the first lifting stage of the transforming element 
xi and x 2 fulfil the equation 



= Pn x 



(81) 



20 



In the following, the first lifting stage 801 is explained, 
which is the lifting stage corresponding to the lifting 

matrix T 3 . 



Let 



be the output vector of the first lifting stage, for 



the time being without rounding to integer values, i.e 
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= 7; 



(82) 



Using the definition of T 3 provided by equation (73), 
equation (82) is re-written as 



2 



I 



AT 
2 



(83) 



Since in this embodiment, a reversible algorithm for integer 
DWT-IV in this provided, there is included a rounding to 
10 integer values. So, according to equation (83), in the first 
step 806 of the first lifting stage 801, xi is multiplied by 
K 3 . The result of this multiplication is rounded to integer 

values in step 807. The rounded values are then added to x 2 
in step 808. Thus, the intermediate signal Zi fulfils the 
15 equation 



(84) 



20 



25 



where [*] denotes rounding operation. 

Since the other four lifting stages 802, 803, 804, 805 of the 
tr ansformin g element illustrated in figure 8, which 
correspond to the matrices T 2 , S, R 2 and R± respectively, 

have identical structures as the first lifting stage 801, 
their explanation is omitted. It only has to be noted that in 
the adding step 809 of the second lifting step 802, Xi 
according to the definition of T 2 , multiplied with D Nn . 
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In the following, the lifting stages of the transforming 
element of the reverse transformation are explained with 
reference to figure 9. 

5 Figure 9 illustrates the lifting stages of the transforming 
element for the reverse transformation of the transformation 
illustrated in figure 8. 

In figure 9, the input of the first lifting stage are the two 
10 blocks of the digital signal y± and z lr z 2 and z$ are 

intermediate signals and xi and £2 are output signals 
corresponding to the first and second block of the digital 
signal, respectively. 

15 The last lifting stage 905 illustrated in figure 9 is inverse 
to the first lifting stage 801 illustrated in figure 8. So, 
in the first step 906 of the last lifting stage 905, Xj is 
multiplied by K s . The result of this multiplication is 

rounded to integer values in step 907. The rounded values are 
20 then subtracted from z x in step 908. Thus, the signal x 2 
fulfils the equation 



Since the other four lifting stages 901, 902, 903, 904 of the 
transforming element illustrated in figure 9, which are 
inverse to the lifting stages 805, 804, 803 and 802 
30 respectively, have identical structures as the last lifting 
stage 905, their explanation is omitted. It only has to be 
noted that after the adding step 909 of the fourth lifting 




(85) 



25 where [*J denotes rounding operation. 



WO 2005/031597 



PCT/SG2004/000122 



39 

step 904 , the result of the adding step 909 is multiplied 
with D Nn to produce Xi- 

It can be seen that the lifting stages 905, 904, 903, 902 and 
901 of figure 9 are inverse to the lifting stages 801 to 805 
of figure 8, respectively. Since also the permutation of the 
input signal corresponding to the matrix can be reversed 

and an according data shuffling stage is comprised by the 
inverse transforming element, the whole provided method is 
reversible, thus, if used in the lossless image encoder 703 
and the lossless image decoder 705 illustrated in figure 7, 
providing a method and an apparatus for lossless image 
coding. 

Although in the explained embodiments, the method according 
to the invention for DCT-IV was used for audio coding and the 
method according to the invention for DWT-IV was used for 
image coding, the method according to the invention for DCT- 
IV can as well be used for image coding, the method according 
to the invention for DWT-IV can as well be used audio coding 
and both can be used as well for the coding of other digital 
signals, such as video signals. 

Considering the equations (63) and (64) one sees that there 
are N/2 roundings in each lifting stage. Therefore, 
considering equation (59), one sees that the total rounding 
nnnihpr i;or the transforming element for the DCT-IV algorithm 
according to the invention is five times N/2, that is 2.5N, 
which is significantly lower than Nlog 2 N according to the 
state of the art. 

Again considering equation (59), one sees that the majority 
of computation power is used in the four N/2 point DCT-IV 
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subroutines corresponding to a multiplication with C™ n , when 

N is a large value, e.g. N = 1024. Because the floating-point 
DCT-IV can be calculated using two half-length DCT-IV plus 
pre- and post-rotations according to the equations (47) and 
5 (4 9) , the arithmetic complexity of the proposed integer DCT- 
IV is roughly estimated to be twice that of the floating- 
point DCT-IV. 

A similar conclusion can be drawn for the integer DWT-IV 
10 transformation function. 



In the following, another embodiment is explained wherein the 
transformation is a discrete Fourier transformation. 



15 The transformation matrix , which is the trans format ion 
matrix of the normalized FFT of order N is given by: 




In = VI l exp x 

m, n = 0, 1, , N - 1 



(86) 



20 where transform size N is even . 



Following the decimation-in-time radix-2 FFT algorithm, In 
can be decomposed in the following way: 



25 



In 



In / 2A/2 In / 
in/ 2 />/2 -in/ 2/V2 



In / 2 



w 



In / 2 



In / 2 



(87) 



Where is, as above, an even-odd matrix, i.e, a 

permutation matrix reordering the components of a vector 
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x = 



x(0) 



x(N - 1) 



by separating the components which correspond to even indices 
from the components which correspond to odd indices, such 



5 that 



x(0) * 
x(2) 

x(N - 2) 
x(l) 
x(3) 

x(N - 1) 



= E 



x(0) 



x(N - 1) 



(88) 



Let In / 2 be the transform matrix of the normalized FFT of 



10 order N/2 . 



Let further W be a diagonal matrix given by 



W = 



4 



< /2 ~ 1 



15 



where 1% = e 



-j2jt / N 



Ijj / 2 denotes, as above, the identity matrix of order N/2 
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In Equation (87), the first matrix from left is an even-odd 
matrix E^o , which only reorders the components in the input 

vector. 



5 In Equation (87) , the second matrix from left can be 

factorized into three lifting matrices in the following way: 



SJ / 2 



IN / 2 



In/ 2 1 [-Q In /2l / 2 

-Q In / 2 IN/2JI 5n/2j|sj/2 In / 2 

(89) 



where Q is an order-JV r /2 permutation matrix given by 



Q = 



0 Jn/2-1 



15 and Jjq / 2-1 i s the order N 12-1 counter, index matrix 



20 



In Equation (87) , the third matrix from left is a counter- 
diagonal matrix, which merely multiplies half of the 
components in the input vector by a complex number residing 
on the unit circle. 



25 



This can be interpreted as a rotation on the complex plane 



Let jc«jc r + jx f be such a component in the input vector 



Let further 



= e -j2kJC / N m cos (2kjt / N) - j sin (2kjc / NJ = 



complex number, i.e. the component of W, which is multiplied 
30 with x when the input vector, after it has been multiplied 
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with the first matrix from the right and the second matrix on 
the right of equation (87) , is multiplied with W. 



The result is 

5 y =* y r + j yi = Wjjjx = (cfcx-r + SkXjJ + j^xi - s k K r) * which equals 

x rotated counter-clockwise by 2kn/N radius on the complex 
plane. Such a rotation can be factorized into three lifting 
steps as follows: 



10 



y r 



-»* c *\[ x '\ 



1 o 

^ i 



1 s 
0 1 



1 0 

SdL o 



(90) 



15 



In Equation (87), the fourth from the right can be factorized 
into three lifting matrices in the following way: 



In / 2/>/2 3n / 
(V2 - 1) In / 2 In / 2 

fir/ 2 



In / 2 - V^2 In / 2 

In / 2 



In / 2 
(V2-i)in/ 2 In/2 



(91) 



The decimation-in-f requency radix-2 FFT algorithm is merely a 
transposition of the decimation-in-time radix-2 algorithm in 
20 Equation (87) . 



Thus, the above procedures also apply to factorising the FFT 
matrix F$ in the decimation-in-f requency way. 



25 Using the factorisation of the right hand side of equation 
(87) in lifting matrices by generating a lifting stage 
corresponding to each lifting matrix a transforming element 
is determined. 
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Since it has been already detailed above how lifting stages 
are generated according to lifting matrices and in this 
embodiment this works analogously to above, the explanation 
5 is omitted here. 
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